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Abstract 

The effect of temporal modulation on traveling waves in the flows in two distinct systems of 
rotating cylinders, both with broken azimuthal symmetry, has been investigated. It is shown 
that by modulating the control parameter at twice the critical frequency one can excite phase- 
locked standing waves and standing-wave-like states which are not allowed when the system 
is rotationally symmetric. We also show how previous theoretical results can be extended to 
handle patterns such as these, that are periodic in two spatial direction. 

PACS numbers: 47.20.-k,47.27.Ak 
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1 Introduction 



The behavior of pattern forming systems under external time-periodic forcing has been the subject 
of intense investigation over the past several years ||^, |2|. Since such periodically forced flows are 
common in nature and technological applications!^, a knowledge of their stability properties may 
have important practical applications. 

There have been several experimental and theoretical studies investigating the effect of 
time-periodic forcing on the stability of the axisymmetric steady spatial patterns produced in the 
flow between two rotating cylinders (Taylor-Couette flow) [Q, ^, ^, |8|, ^, |^]. These studies 
found that temporal modulation of control parameters can stabilize or destabilize the primary 
flows; the spatial properties of the base flow and the axisymmetric Taylor vortex flow remain, 
however, largely unaffected by the modulation. Only a very small shift of the instabilities in the 
parameter space has been observed in these systems. 

For wave structures, however, recent studies by Riecke et al. [jll|] and Walgraef[|l2| 
have shown that for the appropriate frequency of an external forcing a strong resonance may 
occur between the forcing and the spatio-temporal patterns. This can transform initially stable 
traveling waves into standing waves or quasi-periodic structures. This kind of behaviour has 
been experimentally observed in electroconvection of nematic liquid crystals and in binary fluid 



convection by Rehberg et al. [13]. It is also related to the parametric excitation of surface waves 
(Faraday waves) |14[ |. 

In this paper we investigate the effect of temporal modulation on the traveling wave 
patterns that appear near the primary bifurcation boundary of two rotating cylinder systems 
with broken rotational symmetry. First, we will discuss the theoretical background for this study. 
Next, we will describe our experimental studies of the modulated Taylor-Dean and Taylor-Couette 
systems. Finally, we will compare the experimental observations with the theoretical predictions. 

2 Theory 

Using symmetry arguments and suitable amplitude equations it has been shown previously 12, 



15 1 that a resonant temporal forcing can excite standing waves in an effectively one-dimensional 
extended system undergoing a Hopf bifurcation to traveling waves. Here we will extend those 
arguments to systems which are periodic in two directions. This is necessary for spiral vortex flow 
since the azimuthal symmetry of the concentric Taylor system imposes additional constraints on 
the amplitude equations. 

The excitation of standing waves by resonant forcing can be understood intuitively by an 
extension of the standard example of a parametrically forced pendulum. In that case, a periodic 
parametric forcing such as, e.g., the periodic variation of the length of the pendulum, leads to 
a periodic pumping of energy in and out of the pendulum. In order to effect a net increase in 
the pendulum's energy the pumping and the pendulum have to have a certain phase relationship; 
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in particular, the frequencies have to be resonant. Maximal efficiency is achieved if the forcing 
frequency is twice the natural frequency of the pendulum (e.g., decreasing the length of the 
pendulum whenever it travels downward). 

Considering extended systems as an array of coupled pendula makes it clear that only 
standing waves but not traveling waves can be excited parametrically: in a standing wave all 
oscillators have the same phase (albeit different amplitudes), and therefore all of them can satisfy 
the phase relationship necessary for excitation simultaneously; in a traveling wave the phases are 
distributed evenly over the circle and the out-of-phase oscillators will lose the energy gained by 
the in-phase-oscillators. 

To consider spiral vortex flow, which arises in a fluid between counter-rotating concentric 
cylinders |l7|, one has to extend this consideration to two dimensions, since the spirals are periodic 
along the axis and along the azimuthal direction. Therefore, the superposition of oppositely 
traveling spirals ('ribbons') leads to waves which are standing waves along the axis but traveling 
waves in the azimuthal direction. Thus, the periodic forcing may pump energy into the ribbon 
at some azimuthal angle but will extract energy at the same time at a different angle. For a 
net increase these energies have to be different, which is achieved by an azimuthally dependent 
perturbation of the forcing. Alternatively, the azimuthal symmetry of the system itself can be 
broken. Maximal efficiency is achieved, in analogy to the temporal forcing, if the azimuthal wave 
number of the perturbation is twice the azimuthal wave number of the ribbon. 

The above arguments can be made precise by considering spiral vortex flow near its onset 
where its amplitude is small. Then the periodic forcing as well as the azimuthal perturbation can 
also be taken small and the system can be described by coupled amplitude equations. The crucial 
ingredients are the critical eigenvectors of the system, i.e. those modes which have small (or zero) 
growth rate, since they define the center manifold of the system in this parameter regime. In the 
present case the critical modes are the two spiral wave modes of the flow field v as well as the 
periodic forcing F and the azimuthal perturbation P, 

V = A{T) e^''^+*'^h*+^'^<^ f„(r) + B{T) e^i'-^'^^t-imcp f^^^^ ^ ^ ^ ^ ^^^^^^^ 

F = be''^^^ + c.c, (2) 

p = we'^'^ + c.c. (3) 

Here T is a slow time and ^a,b{T~) are the radial eigenvectors of the spiral vortex flow. The axis 
of the cylinders is taken along the z-direction and (/> is the azimuthal angle. The frequencies LVh 
and Ue are the Hopf frequency of the spirals and the frequency of the forcing, respectively. Below 
we will choose u)e close to 2liJ/j. The difference will determine the detuning. Although the forcing 
and the azimuthal perturbation are imposed externally we consider them as dynamical variables 
and include them in the center manifold. This simplifies the derivation of the relevant amplitude 
equations for A{T) and B{T). Their form is obtained by considering the most general polynomial 
which is consistent with the symmetries of the system (in the absence of forcing, F = 0, and the 
azimuthal perturbation P = 0): translation symmetry in space along z, rotational symmetry in 
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(/), translational symmetry in time t, and reflection symmetry in z. These symmetries allow a 
linear coupling of the left and right spiral amplitudes via the forcing and azimuthal perturbation 
only if the frequencies and ujfi as well as the azimuthal wave numbers m and n are related as 

^ = k, ^ = 1, k,l = 1,2,3... (4) 

UJe n 

Simply put, under these conditions the term B ^iqz-iu}ht-im<t>^ikujet^iin<f) corresponds to the 
same Fourier coefficient in z, t and (p as ^g292+«t^ht+«m,0 ^^^^ ^^^^ therefore arise in the evolution 
equation for A, thus providing the crucial linear coupling. To cubic order one then obtains the 
following equations, 

drA = aA + d6V™/"S + c|^|2^ + 5|SpA, (5) 
drB = a*B + db''w^"'/''A + c*\B\'^B + g*\A\'^B, (6) 
= 0, dTW = 0. (7) 

The real part of the coefficient a gives the growth rate of the spiral waves and is proportional 
to the distance from the Hopf bifurcation, a^ = a{Ri — Ric), with Ri the Reynolds number for 
the inner cylinder and Ric its critical value. To allow for a small detuning in the forcing, the 
frequency Uh in ( |I|) is replaced by a;e/2. The detuning ^ is then related to the imaginary part Oj 
as 

CJft - UJe/2 tti - 70^ , , 

? = = , (oj 

where the term 7a,. accounts for the change of the linear frequency with the control parameter 
Ri. The coefficients d, c and g are in general complex whereas the product db^uP'^/^ can be 
chosen real. All coefficients are functions of the invariants |6p and This is most relevant 

for the small coefficient a. Thus the forcing and the azimuthal perturbation lead to a shift in the 
threshold and of the Hopf frequency even in the non-resonant case. 

Eqs.( |5|, ^) show that the forcing and the azimuthal perturbation have maximal effect 
if /c = 1 and / = 2m /n = 1, i.e. if the external frequency is twice the Hopf frequency and the 
wave number of the azimuthal perturbation is twice that of the spiral waves. Unfortunately, 
the latter condition cannot be implemented easily, since eccentrically mounted cylinders yield a 
perturbation of the form cos(/>, i.e. n = 1. Clearly, in the absence of any azimuthal perturbations 
(w = 0) no linear coupling of oppositely traveling waves occurs and the standing waves in question 
cannot arise directly from the basic state A = Q = B through a linear instability. 

The amplitude equations ( ^, ^) have been previously analyzed in detail jl^, 12 1 for the 



case gr < Cr < 0, i.e. for a Hopf bifurcation to stable traveling waves rather than (unlocked) 
standing waves. The main result is the excitation of standing waves, \A\ = \B\ = const, by 
the periodic forcing below the Hopf bifurcation (a^ < 0). As is clear from the fact that the 
amplitudes A and B are constant in time these waves are phase-locked to the forcing (cf. (|l]) 
with uJh — > uJe/2). A sketch of a typical phase diagram is shown in Fig. 1 where the loci of various 
transitions are given as functions of the forcing strength (measured by b) and the mean Reynolds 
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number (measured by ar). The transition from the basic state of purely azimuthal flow to the 
phase-locked standing waves (ribbons) occurs along the line marked PSW, which is given by 



Along the line marked H the basic state undergoes a Hopf bifurcation to traveling waves as well as 
(unstable) standing waves. This bifurcation exists already in the absence of any periodic forcing 
(5 = 0). Along the line marked PB the phase-locked standing waves become unstable to the 
traveling waves in a secondary parity-breaking bifurcation. Not relevant for this experiment are 
the transition of the unstable phase-locked standing waves to unstable standing waves at SW and 
the saddle-node bifurcation of the phase-locked standing waves along SN. 

To obtain quantitative results pertaining to spiral vortex flow the linear as well as the 
nonlinear coefficients have to be determined numerically. This is a formidable task and will not 
be attempted here. It should be noted, however, that to leading order the nonlinear coefficients 
c and g are not affected by the forcing and the azimuthal perturbation. Thus, the values for the 



concentric case can be used, which have been determined previously [16|. 

In addition to the results on modulated spirals in the slightly eccentric Taylor system 
we will also present results on traveling inclined rolls, which occur in the Taylor-Dean system [ 24 1 



In this system the azimuthal symmetry is broken strongly. Therefore, instead of the discussion 
presented above the original derivation of the amplitude equations in the presence of a single 
translation symmetry applies [11, 12, 15 1. It is based on an expansion of the flow field in the form 

V = A(r) 6*9^+^'^'^* fa(r, 4>) + B(T) e*?^-^'^^* %(r, (j)) + c.c. + h.o.t. (10) 

and leads to (^J6|) with w = I. 



3 Experimental Apparatus and Procedures 

Two variations of the concentric rotating cylinder (Taylor-Couette) system were employed to test 
the concepts of section 2. The systems break the rotational symmetry in very different ways. 
In the next few sections, we will describe each one in detail. First, we present the experimental 
system with strongly broken rotational symmetry, the Taylor-Dean system. 



3.1 Taylor-Dean System 

The Taylor-Dean system consists of two independently rotating horizontal coaxial cylinders with 
a partially ffiled gap (see Figure 2)p3|, p^ j. The partial filling of the gap in the Taylor-Dean 
system breaks the rotational symmetry of the flow. The rotation of the cylinders and the two free 
surfaces impose a pressure gradient along the azimuthal direction. As a result, the flow sufficiently 
far away from the free surfaces is a combination of Couette flow due to rotation of the cylinders 
and Poiseuille flow due to the azimuthal pressure gradient. The main control parameters of this 
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system are the inner and outer cylinder Reynolds numbers, Ri = 2Tifirid/u and Rq = 2ti foVod/v, 
where fo and fi are the outer and inner cylinders' rotational frequencies, d is the gap width, 
and v is the kinematic viscosity of the fluid. As the control parameters are varied the base flow 
instabilities change from those associated with Taylor-Couette to those associated with Dean flow. 
The first experimental work with the Taylor-Dean system was done by Brewster and Nissan [p2| , 
who studied the threshold of the instability when the inner cylinder rotates and the outer cylinder 
is stationary. More recently, Mutabazi, Normand, Peerhossaini,and Wesfreid p5[| have solved the 
linear stability problem for axisymmetric and non-axisymmetric perturbations in the flow. They 
found both stationary and traveling wave instabilities depending on the ratio of the angular 



velocities /i (= fo/fi )|23|. We kept the outer cylinder stationary {Rq = 0) for this study, while 
rotating the inner cylinder. In this case the flrst transition from the unperturbed base flow is to 
traveling inclined rolls as the inner cylinder speed increases, as shown in the experimental study 
of Mutabazi, Hegseth, Andereck and Wesfreid[^]. The transition to traveling inclined rolls is, 
within the experimental error, a supercritical Hopf bifurcation. At onset these rolls also have 
no preferred direction and may move either left or right along the cylinder axis. Therefore, this 
system can be used to verify the theoretical prediction that a breaking of the time translational 
symmetry by a small periodic modulation of the control parameter will result in a stable standing 
wave pattern[0, p^ . 

Our experimental system consists of an inner cylinder made of black Delrin plastic with 
radius rj = 4.47cm and a stationary outer cylinder made of Duran glass with radius Tq = 5.08cm, 
giving a gap d = — ri = 0.59cm and radius ratio r] = ri/rg = 0.883. Two plastic rings are 
attached to the inner cylinder a distance L = 52.4cm apart, giving an aspect ratio T = L/d = 88, 
large enough to conceive of this as an extended system where one can neglect end effects. In this 
system the fllling level fraction n = 6f/2TT has been fixed at 0.75, where Of is the filling angle. 

The working fluid was pure double distilled water or a solution of double distilled water 
and 44% glycerol by weight with 1% of Kalliroscope AQIOOO added for visualization. These 
nearly two-dimensional (~ 30 x 6 x 0.07^m) polymeric flakes align along the streamline surfaces, 
reflecting light according to their orientation. Generally the dark areas indicate flow along the 
observer's line of sight, while the light areas indicate flow perpendicular to the line of sight. The 
apparatus was kept in a temperature controlled room so that the temperature of the working 
fluid was held constant to within O.l'^C. 

Spatial and temporal properties of the flow patterns were obtained using a 512 x 480 
pixel CCD camera connected to an image processor. The image processor board, installed in a 
PC, captured a picture of the flow pattern and then a software routine was used to obtain the 
intensity along a single line parallel to the axis of the cylinders. This system is able to process up 
to one line every 0.11 sec. The data were then transferred to a VAX 4000-90 computer system and 
an analysis of the resulting intensity versus axial position as a function of time plots (space-time 
data plots) yields the wavelengths and the dynamics of the patterns in time and space. 

In the Taylor-Dean system, we have used a combination of two stepper motors (Com- 
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pumotor A83-93) to drive the inner cylinder, one motor to produce a net rotation and the other 
to produce a sinusoidal modulation of the inner cylinder angular velocity. The first motor was 
directly connected to the inner cylinder and rotated with constant angular velocity. The housing 
of this motor was oscillated by a push-rod arrangement mechanically coupled to the second mo- 
tor (see Figure 2). The rotation of the second motor gave a net output at the inner cylinder of a 
constant angular velocity plus a periodic sinusoidal variation in angular velocity. 

The motor speeds were controlled through Compumotor 2100 Series indexers and could 
be changed either manually or through computer control. The motor speeds have a frequency 
accuracy of 0.02%. Typically, a computer controlled the rotation speed, direction, and ramp- 
ing rates. Because both the frequency and amplitude of the modulation could be varied this 
introduced two new control parameters into the system, as required by the theory. 

The Reynolds number of the inner cylinder is now R = Ri + RmSin{2TTfmt), where Rm = 
iTrfaTid/v and fa is the amplitude of the sinusoidal rotational frequency. The two dimensionless 
parameters are Rm/Ric and the detuning parameter ^ = {2fh — fm)/'^fh, where Ric is the critical 
Reynolds number for onset of the traveling roll pattern and fh is the frequency of the teaveling 
rolls (Hopf frequency) in the absence of modulation. 

To obtain the location of the onset of patterns for each modulation frequency fm and 
modulation amplitude Rm, we employed the following method. We first set the Ri value below 
the onset of instabilities and then increased it quasistatically (keeping both amplitude Rm and 
frequency fm of modulation at a fixed value) until a flow pattern appeared in the system. Then 
several sets of space-time data were taken while increasing Ri. We repeated this procedure with 
increasing amplitudes and various frequencies around twice the Hopf frequency (the frequency of 
the traveling rolls near onset when there is no modulation). 

3.2 Results and Discussion 

When there is no modulation {Rm = 0), the base flow bifurcates supercritically to a traveling roll 
pattern as we increase the inner cylinder speed. A typical space-time diagram for the traveling 
roll pattern near onset (i?j = 265, Rq = 0, and Rm = 0) is shown in Figure 3. The wavelength of 
the rolls along the cylinders is A = 0.841 cm. At e (= ^'j^.^'" ) slightly greater than 0.0 the pattern 
fills most of the working space and both left and right traveling rolls may exist with a vertical 
defect line between them (see Figure 4(a)). The frequency of the traveling rolls (Hopf frequency) 
near the onset ( Ric = 263) is 0.543 Hz. Upon further increase of Ri, the flow undergoes a second 
instability to a short wavelength modulation of the traveling rolls at Ri = 303, and then to an 
incoherent pattern at about Ri = 338 [p7|]. 

The transition sequence changes dramatically with modulation of the inner cylinder 
speed. When we modulated the inner cylinder sinusoidally near detuning parameter ^ ~ 0, we 
found standing waves rather than the traveling rolls (see Figure 4(b)), as predicted by theory. 
Figure 5 shows a phase diagram of the primary transitions from the base flow to standing waves 
and the secondary transition from standing waves to traveling waves for ~ as we varied 
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the modulation amplitude Rm and inner cylinder Reynolds number One interesting feature 
to note in Figure 5 is that, when the modulation amplitude is increased above a critical value 
{Rm/Ri > 0.05), the standing waves can be excited at Ri values much lower than the critical Ri. 
Another interesting feature is that as the amplitude of modulation increases, we observe standing 
waves over a widening range of Ri. The traveling wave state reappears when Ri increases, as has 
been observed in other systems [p!^] and in agreement with the theory 

At small modulation amplitudes {Rm/Ri < 0.05) only the traveling roll state appears 
when Ri is increased, also in agreement with the theoretical predictions (c.f. Figure 1). 

The space-time diagram of the standing waves for Rm/ Ric = 0.30, Ri = 241 and ^ = 
lb 0.01 is shown in Figure 6. During one modulation period the light intensity at a given axial 
position varies periodically, indicating the presence of standing wave patterns. The frequency of 
the standing waves is half the modulation frequency. 

A quantitative analysis of the patterns was carried out using 2-dimensional Fourier 
transforms of the patterns in time and space. This proved to be a very useful method for 
decomposing the standing wave patterns into their left and right traveling wave components from 
the original space-time CCD data. From the decomposition we obtained the spatial wavelengths, 
temporal frequencies, and the amplitudes of each component. Figures 7(a) and 7(b) shows the 
right and left traveling waves obtained from the space-time data shown in Figure 6. The power 
spectra obtained using the 2-dimensional FFT for both right and left traveling waves are shown in 
Figures 7(c), and 7(d). The left and right components have similar space and time characteristics. 
The small differences in amplitudes of the spectral peaks may be attributed to slightly nonuniform 
lighting conditions. 

At a still higher inner cylinder speed the standing waves lose their stability to a trav- 
eling wave state. This transition is supercritical within our experimental resolution. No mixed 
standing/traveling rolls states have been seen. The space-time diagram for the traveling wave 
state just above the onset at Ri = 269, Rm/Ric = 0.212 is shown in Figure 8. The decomposition 
of this space-time data and associated power spectra are shown in Figures 9(a), (b), (c), and (d). 
These figures show that one traveling mode dominates over the other, much weaker, components 
present in the system. 

We also tested the sensitivity of the onset of the patterns, and which pattern appears 
first, to variation of the detuning parameter ^ in this system. The result for Rm/Ri = 0.3, shown 
in Figure 10, indicates that the onset of the patterns is very sensitive to the value of the detuning 
parameter ^ for this system. When ~ the standing wave pattern appears at Ri well below 
Ric = 263. As we changed ^ away from ^ = (lower or higher modulation frequencies than 2//j) 
the onset of the patterns was delayed considerably. In fact, at higher modulation frequencies 
< —0.2) patterns appeared only at supercritical Ri{> Ric) values for Rm/Ric = 0.3. Figure 10 
also shows that, for fixed Rm, as we move away from ^ = the characteristics of the patterns also 
change, so that traveling wave patterns appear at onset (^ < —0.08 or ^ > 0.04 for Rm/Ric = 0.3) 
rather than standing waves. It is possible that even for these large values of |4| standing waves 
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would reappear as the critical mode for higher fixed Rm/Ric than 0.3. 

Figure 11 shows the ratio of the fractional power under the left (the one induced by 
the modulation in this case) traveling wave spectra to the total power, as a function of ^. This 
characteristic resonance curve illustrates the disappearance of the left traveling wave as the de- 
tuning parameter ^ shifts away from 0. The maximum fractional power was a little less than 
the expected 0.5, which can be attributed to the non-ideal lighting conditions. We also analysed 
the frequencies of the two components as a function of ^. Figure 12 shows the frequency ratio 
as a function of This plot shows the standing wave region clearly, where the two components 
lock to the same frequency over a range of ^. The primary transition to standing waves and 
the secondary transition to traveling waves in Rm / Ric Ri parameter space are very sensitive 
to ^ for ^ near zero. Figures 13 (a) and (b) give the phase diagrams when ^ is equal to —0.03 
and -1-0.03, which show these phase changes clearly. A comparison of Fig. 13 and Fig. 5 shows 
that the minimal modulation amplitude for the excitation of standing waves is actually lower for 
^ = 0.03 than for ^ = 0. This can be ascribed to the dependence of the linear frequency of the 
waves on Ri as well as j6p and \w\^ as expressed in (^; optimal resonance is obtained for = 
rather than ^ = (cf. ®). 

Further increases in Ri beyond the traveling wave state resulted in the appearance of 
short wavelength modulations near Ri ~ 300 and incoherent patterns at Ri ~ 340 for Rm/Ric = 
0.25 and ^ = 0. These patterns and their decomposed left and right traveling components are 
shown in Figures 14 and 15. 



3.3 Eccentric Taylor-Couette System 

In the Taylor-Dean system, the rotational symmetry is strongly broken by the air/liquid interfaces. 
We now consider a system with slightly broken rotational symmetry, which can be modeled as a 
perturbation on a system with two translational symmetries (strictly speaking, a pure translation 
in space along z and a rotational symmetry along 0). 

The classical concentric Taylor-Couette system has translational symmetry in the axial 
direction and rotational symmetry around the axis. As we discussed in section 2, for this system 
it is necessary to apply periodic forcing as well as an azimuthal perturbation to linearly couple 
the left and right traveling spirals to produce the standing waves. We produced an azimuthal 
perturbation in this system by offsetting the axis of the inner cylinder relative to that of the outer 
cylinder while maintaining the two axes parallel, thereby making the system eccentric. 

The experiment was performed in the region of (i?o, Ri) parameter space where traveling 
waves in the form of spiral vortices occur as the primary instability. This happens when the 



cylinders are counter-rotating and concentric Jl7| , [Lq , 19, The spiral patterns travel in both 
the axial and azimuthal directions, and they break both the axial and azimuthal symmetry of the 
base flow. For a given radius ratio there is a unique value of outer cylinder speed above which 
the primary bifurcation from the base flow is a supercritical Hopf bifurcation to the time periodic 



spiral flow|15, 19|. The azimuthal wave number (m) of the spirals increases as the outer cylinder 
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speed increases. The locations of the crossover points between spirals with different azimuthal 
wave numbers are uniquely determined by the radius ratio of the two cylinders [ 19 1. 

Our eccentric Taylor-Couette cylinder system, shown schematically in Figure 16, consists 
of an inner cylinder made of black Delrin plastic with radius rj = 4.76cm, and an outer cylinder 
made of Plexiglas with radius Tq = 5.95cm, which gives a gap d = — ri = 1.19cm, and a radius 
ratio T] = 0.800. The main control parameters are the inner and outer cylinder Reynolds numbers, 
Ri = 2-Kfirid/v and Ro = 2t: foVod/v, where d is the average gap width. A stationary support 
holds the apparatus at the top and the bottom. The outer cylinder is connected to the stationary 
support by means of bearings at both ends. The upper end of the inner cylinder is attached to a 
long shaft which hangs from a horizontally movable plate on the stationary framework (see Figure 
16). The lower end of the inner cylinder is left unattached, and is suspended « 1mm above the 
bottom of the system. Eccentricity is adjusted by offsetting the axis of the inner cylinder relative 
to the fixed axis of the outer cylinder. The position of the axis of the inner cylinder is read to 
an accuracy of 0.005cm using a micrometer attached to the stationary support. The eccentricity 
is e = £/{to — ri), where e is the offset of the two cylinder axes. To maintain consistent end 
conditions, the upper and lower boundaries of the flow are formed by Teflon rings attached to 
the outer cylinder and located near the ends of the cylinder. There is a narrow gap of 0.4 cm 
between each of these rings and the inner cylinder that is nevertheless wide enough to allow for 
offsetting the inner cylinder. The length of the fluid column L is 40.40 cm, giving an aspect 
ratio r = = 34.0, large enough to minimize end effects. Again, the working fluid was 

pure double distilled water or a solution of double distilled water and 44% glycerol by weight 
with 1% of Kalliroscope AQIOOO added for visualization. The experimentally obtained critical 
Reynolds number for concentric cylinders when = has been compared with reference values 
to check the accuracy of u obtained from tabulated data ||2l[| . The difference was less than 1%. 
Both inner and outer cylinders are driven by two independently rotating Compumotor stepper 
motors, which are controlled through Compumotor 2100 Series indexers. Since the modulation 
amplitude required in the eccentric Taylor-Couette system was much greater than for the Taylor- 
Dean system, we could not use the two motor push rod arrangement to produce both linear and 
sinusoidal modulation of inner cylinder. Therefore, we have used a single motor to produce both a 
constant rotation speed component and the sinusoidal modulation of the inner cylinder rotation 
rate by sending control commands (a sequential change of shaft motion parameters with time 
delays) to the indexer using the computer. This produces an approximate sinusoidal variation, 
with more than 40 step changes in velocity per cycle. 

To obtain the location of the onset of patterns for each modulation frequency fm and 
modulation amplitude Rm, we employed a method similar to that used in the experiment on 
the Taylor-Dean system. We kept the outer cylinder speed fixed at a value within the range 
where the bifurcation to spiral vortices occurs. Then, at a fixed modulation amplitude Rm and 
constant frequency fm, we increased the average inner cylinder speed until the pattern appeared 
for fixed modulation amplitude and constant frequency. Then several sets of space-time data 
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were taken while increasing Ri. We repeated this procedure with increasing amphtudes Rm 
and various detuning parameters ^ (i.e., frequencies around twice the Hopf frequency, where the 
Hopf frequency in this case is the frequency of the spiral patterns near onset in the absence of 
modulation) . 



3.4 Results and Discussion 

When there is no modulation the base flow bifurcates supercritically to spirals as we increase 
the inner cylinder speed while keeping the outer cylinder speed constant. A typical space-time 
diagram for the spiral pattern near onset {Ri = 151, Rq = 155, and Rm = 0) is shown in Figure 
17. Here the intensity of a line parallel with the axis of the cylinders was recorded every 0.6s for 
100s. The wavelength of the spirals along the axis of the cylinder is A = 2.07cm ~ 2d. Light sheet 
visualization through the gap cross section shows that the spirals exist near the inner cylinder. At 
g (= ^i-^ic ) slightly greater than zero the pattern fills most of the working space along the axis. 
An upward moving spiral exists near the bottom of the system and a downward moving spiral 
exists near the top. A horizontal defect line forms where the two spirals meet. These defect lines 
are not necessarily halfway between the top and the bottom along the axis of the cylinders. Such 
defects are inherent to traveling wave patterns [ 26 1. The frequency of the spirals (Hopf frequency) 
with azimuthal wavenumber m = 2, as shown in Figure 17, near the onset {Ric = 151, i2o = 155 
for e = 0) is O.lUHz. 

When we sinusoidally modulated the inner cylinder rotation speed, we observed wave 
patterns which resemble standing waves (time-dependent patterns which are stationary in space) 
rather than traveling waves (spirals) as the first bifurcation from the base flow. Figure 18 shows a 
space-time diagram of the standing-wave-like pattern near the transition to the spiral state with 
azimuthal wave number m = 2 at i?j = 151, Rq = 155, ^ = 0.4, e = 0.126, and ^ 0. Figure 
19(a) and 19(b) shows the resulting right and left traveling waves obtained from the space-time 
data shown in Figure 18. The power spectra obtained using the 2-dimensional FFT for right 
and left traveling wave components, shown in Figures 19(c) and 19(d), have the same frequency 
characteristics, but different peak amplitudes. A true standing wave would have equal amplitudes 
in the traveling components. Even at very high modulation amplitudes Rm, one spiral component 
always dominates the other. Figure 20 shows the ratio of the amplitudes of the primary peaks of 
the two spiral components versus modulation frequency fm- 

Figure 21 shows the phase diagram of Rm/Ric vs Ri for e = 0.168 and ~ 0. This shows 
that the base flow became unstable at lower Ri as the amplitude of modulation Rm increased, in 
agreement with the theoretical predictions. Also, the standing-wave-like patterns appeared only 
at large modulation amplitudes ( for e = 0.168, > 0.17) when the detuning parameter ~ 0. 
At small modulation amplitudes {Rm/R-i < 0.17) only the traveling roll (spiral) state appears 
when Ri is increased. Whether standing- wave-like patterns are produced also depends on the 
eccentricity of the system. For small eccentricity e, they appeared only at higher modulation 
amplitudes, in agreement with the theoretical predictions. 
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Compared with other systems (e.g., the Taylor-Dean system discussed in the previous 
section or electroconvection of nematic liquid crystals and binary fluid convection [l^), in this 
system, the standing-wave-like patterns with two components having the same frequency appeared 
only in a small region of parameter space near the resonance frequency (see Figure 22). Both 
the small region of parameter space and the unequal power spectra amplitudes indicate that the 
coupling of the modulation to the fluid flow is very weak. Two factors may contribute to this. 
Eqs.( |5|, ^) show that the azimuthal symmetry breaking employed in the experiment (n = 1) 
enters the coupling between left- and right-traveling spirals only with its fourth power for spirals 
with m = 2. Thus, quite large eccentricities are needed to give a strong effect. Motivated by 
this observation we have also studied the modulation of m = 1-spirals. However, no resonant 
excitation of standing waves was observed; instead, the patterns relax to the base flow during 
part of the cycle. Since the frequency of the waves is very small, they appear to follow the 
periodic forcing adiabatically indicating that their growth rate is of the same order as the forcing 
frequency. Therefore for ( ^, ^) to apply, the growth rate would have to be reduced substantially. 
This would require much higher resolution in the average Reynolds numbers than is accessible in 
the present apparatus. The investigation of m = 1-spirals is further complicated by the presence 
of other modes like axisymmetric vortices and interpenetrating spirals at nearby parameter values. 

A second reason for the apparent weakness of the coupling is presumably the penetration 
depth of the oscillation into the bulk of the fluid, as characterized by the viscous Stokes layer 



of width 6 = y/u/uj, near the inner cylinder. The thickness of the Stoke's layer 5 at ^ = for 
both m = 1 and 2 spiral states are 0.24cm and 0.17cm, much smaller than the size of the gap 
d = 1.19cm. Thus, the bulk of the fluid feels the modulation only very weakly. (It is possible that 
similar effects may have been present in the binary fluid convection experiment reported in [13|. 
This may have contributed to the less well-defined standing waves reported there, in contrast 
with the electroconvection experiments.) Futhermore, in the eccentric Taylor-Couette system, 
modulation creates only a weak azimuthal pressure gradient due to the small eccentricity, and 
the oscillatory Stokes layer influence may not be enough to excite a strong secondary traveling 
component. In the Taylor-Dean system, on the other hand, the inner cylinder modulation cre- 
ates not only the oscillatory boundary layer influence but also an induced oscillatory azimuthal 
pressure gradient across the whole gap, due to the free surface, which more readily affects the 
pattern. 

The standing-wave-like patterns appearing in the eccentric Taylor-Couette system at 
the ^ = near the transition boundary to m = 2 spiral state lose their stability to an incoherent 
pattern when the inner cylinder speed is increased further above the onset. This may be due 
to the loss of stability of the basic spirals to interpenetrating spirals and wavy spirals. Even for 
the concentric cylinder case (e = 0) we have seen some evidence of a weak secondary spiral at 
very high modulation amplitudes. This could be due to small imperfections in the experimental 
apparatus (e.g. nonparallel cylinder axes) which may break the rotational symmetry sufficiently 
to induce a second spiral. 
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4 Conclusions 



In conclusion, we have found that time-periodic modulation at close to twice the frequency of 
a Hopf bifurcation can induce standing waves in systems with traveling wave patterns in two 
directions. The Taylor-Couette system with counter-rotating cylinders produces spirals at onset 
over a large parameter range, but unless the rotational symmetry of the apparatus is broken 
standing waves cannot be excited by modulation. This is in contrast with the case of one- 
dimensional traveling wave patterns studied in convection [p!^], where only temporal modulation 
is necessary to excite standing wave patterns. Temporal modulation of the eccentric Taylor- 
Couette system induces a second traveling spiral pattern as predicted by theory. However, owing 
partly to weak coupling of the oscillations to the bulk of the flow, it never grows to a large enough 
amplitude, nor does it couple with the other traveling spiral sufficiently strongly, to produce simple 
standing wave patterns at onset. The temporally modulated Taylor-Dean system, in contrast, 
with its strongly broken azimuthal symmetry, yields clear agreement with the qualitative features 
of the theoretical model. 
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Figure Captions 

Figure 1: A sketch of a typical phase diagram for modulated rotation at resonance (uje = 2ujh), 
where the loci of various transitions are given as a function of the forcing strength (measured by 
6) and mean Reynolds number (measured by a^). The transition from the basic state to standing 
waves occurs along the line marked PSW and from the basic state to traveling waves occurs along 
the line marked H. 

Figure 2: (a) Schematic diagram of the Taylor-Dean apparatus, (b) Schematic cross section of the 
apparatus. The front side is the side where the observer sees the inner cylinder rotating upward 
as shown. 

Figure 3: Space-time diagram of the traveling roll pattern in the Taylor-Dean system when 
Ri = 265, Ro = 0, and Rm = 0. The wavelength of the pattern is A = .841cm and frequency = 
0.543 Hz. 

Figure 4: Pictures of the (a) traveling roll (for Ri = 269 and R„i = 0) and (b) standing wave 
states in the Taylor-Dcan system (for Ri = 265 and Rm/Ric = 0.28). 

Figure 5: Phase diagram of standing wave and traveling wave states of the Taylor-Dean sys- 
tem when the detuning parameter is ^ = 0. The plus and cross signs indicate the location of 
the bifurcation from the base state to standing waves and standing waves to traveling waves 
respectively. 

Figure 6: Space-time diagram of the standing wave state at Rm/Ri = 0.30, Ri = 241 and 
^ = ± 0.01 (Taylor-Dean system). 

Figure 7: Space-time diagrams of (a) the right traveling wave component and (b) the left traveling 
wave component of the standing wave state shown in Figure 6, obtained using a 2-dimensional 
FFT decomposition, (c) and (d) are the frequency power spectra of the decomposed right and 
left traveling modes. The frequencies are scaled to the inverse of the diffusive time scale u/d^. 
Figure 8: Space-time diagram of the traveling wave state at Rm/Ric = 0.212, ^ = 0.0, and 
Ri = 269 just above the onset (Taylor-Dean system). 

Figure 9: Space-time diagrams of (a) the right traveling wave component and (b) the left traveling 
wave component of the traveling wave state shown in Figure 8. (c) and (d) are the frequency 
power spectra of the decomposed right and left traveling modes. The frequencies are scaled to 
the inverse of the diffusive time scale z^/d^. The right traveling wave component dominates the 
left traveling component in the system. 

Figure 10: The onset of the primary flow transition Ri vs the detuning parameter ^ (= {2fh — 
fm)/'2^fh) at Rm/Ric = 0.3 (Taylor-Dean system). 

Figure 11: The fractional power of the secondary traveling wave state vs the detuning parameter 
^ at Rm/ Ric = 0.3. 

Figure 12: The frequency ratio of the two traveling wave components vs the detuning parameter 
^ at Rm/Ric = 0.3. Both left and right traveling components have the same frequency within 
experimental error over a fairly wide region near resonance ^ = 0. 
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Figure 13: Phase diagrams when (a) $, = —0.03 and (b) ^ = 0.03. Pluses and crosses indicate the 

transitions from the base flow to standing waves and to travehng waves respectively. 

Figure 14: (a) The traveling wave state with the short wavelength modulations at Ri = 300 and 

Rm/Ric = 0.254 (Taylor-Dean system), (b) and (c) arc the decomposed space-time data for right 

and left traveling modes, respectively, of the same traveling wave state shown in (a). 

Figure 15: (a) The incoherent patterns at Ri = 340 and Rm/Ric = 0.254 (Taylor-Dean system). 

(b) and (c) are the decomposed space-time data for right and left traveling modes, respectively, 

of the incoherent pattern shown in (a). 

Figure 16: (a) Vertical, and (b) horizontal, cross-sections of the eccentric cylinders system. 
Figure 17: Space-time diagram of the spiral pattern (m = 2) in the Taylor- Couette system at 
Ri = 151, Ro = 155, e = 0, and Rm = 0. 

Figure 18: Space-time diagram of the standing-wave-like pattern in the Taylor-Couette system 
at Ri = 151, Ro = 155, e = .126, Rm/Ric = 0.4, and ^ = 0. 

Figure 19: (a) and (b) arc the decomposed components of the space-time data in Figure 18. (c) 
and (d) are the frequency power spectra of the right and left traveling modes. The frequencies 
are scaled to the inverse of the diffusive time scale u/d"^. The right traveling wave component 
dominates the left traveling component in the system. 

Figure 20: The ratio of the spectral peak amplitude of the secondary spiral wave state to the 
spectral peak amplitude of the primary spiral wave state vs the detuning parameter ^ at i?i = 
151, Ro = 155, Rm/Ric = 0.53, and e = 0.126. 

Figure 21: Phase diagram of Rm/Ric versus Ri for the eccentric Taylor-Couette system with 
eccentricity e = 0.168 and detuning parameter ^ = 0. 

Figure 22: The frequency ratio of the two traveling wave components versus the detuning pa- 
rameter in the eccentric Taylor-Couette system at Ri = 151, Rq = 155, Rm/Ric = 0.53, and 
e = 0.126. The left and right traveling components have the same frequency only in a very small 
region near resonance ^ = 0. 
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